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The Rb-NH interaction is investigated as a prototype for interactions between alkali metal atoms 
and stable molecules. For spin-aligned Rb and NH that interact on a quartet surface {^A"), the 
interaction is relatively weak, with a well depth of 0.078 eV. However, there are also doublet surfaces 
of ion-pair character that are very much deeper (well depth 1.372 eV) and may be important for 
atom - molecule collision rates. Similar deeply bound ion-pair states are likely to exist for other 
alkali atom - molecule pairs. It is shown that a transition to a dipolar superfluid phase will be easier 
to achieve for low-mass than for high-mass species. 
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There is great interest in the production and trap- 
ping of translationally cold molecules, and in the pos- 
sibility of achieving quantum degeneracy in molecular 
gases P, 12 0j B- There is special interest in dipolar 
quantum gases [5j, where the trap geometry controls the 
attractive or repulsive nature of the anisotropic dipole- 
dipole interaction. Important goals in this area include 
Bose-Einstein condensation of dipolar gases 0], realiza- 
tion of phase transitions of such gases in optical lattices 
0, and formation of superfluid pairs and observation of 
the Bardeen-Cooper-Schrieffer (BCS) transition to a su- 
perfluid state in dipolar Fermi gases The electric 
dipole moments of ultracold polar molecules could also 
be used as qubits in quantum computation 0. 

One very promising route for the production of cold 
dipolar molecules is molecular beam deceleration [lfij |. 
Meijer and coworkers have developed a decelerator based 
on the molecular Stark effect and switched inhomoge- 
neous electric fields and have used it to slow a 

beam of ND3 molecules and trap them at a temperature 
around 25 mK [11. Beam deceleration is also applicable 
to other polar molecules, such as OH and II2O, but is 
typically only likely to achieve temperatures of the order 
of 1 niK. To cool the molecules further, towards tem- 
peratures where quantum degeneracy might be achieved, 
another cooling technique is needed. One possible route 
is sympathetic cooling b y t hermal contact with a cold 
atomic gas such as ^^Rb 

The purpose of the present paper is to investigate 
the potential energy surfaces for interaction of a polar 
molecule with an alkali atom, in order to help under- 
stand collisions in molecule - alkali atom mixtures. For 
this prototype study, we choose NH interacting 
with with Rb i^S). This system is topical because Mei- 
jer and coworkers ;14j have recently proposed a scheme 
for the production of ultracold NH molecules in their 
ground electronic and vibrational state. In this scheme, 
NH molecules in their long-lived a^/S. metastable state 
are first slowed in a Stark decelerator. The molecules are 
then excited to the w = level of the A^Ti excited state. 



allowed to decay spontaneously to the v — level of the 
X^ll~ ground state, and finally trapped magnetically. 

An important and general point is that many molecule 
- alkali atom systems have deeply bound electronic states 
with ion-pair character. Alkali metal atoms have low ion- 
ization energies, while molecules containing electronega- 
tive elements such as C, N, O and the halogens can ac- 
cept an extra electron relatively easily. There is thus 
an ion-pair state (Rb"*" NH~ in the present case) that is 
asymptotically not far above the energy of the neutrals. 
Because of the Coulomb attraction between the ions, the 
potential surface for the ion-pair state cuts down fast in 
energy with decreasing separation, and around the equi- 
librium geometry is often the lowest electronic state. As 
will be seen below, this is the case for RbNH. 

The X^'Ti~ ground state of NH has 2 electrons with 
parallel spins in tt antibonding orbitals. At coUinear ge- 
ometries (with Coov point group symmetry), there are 
two electronic states of RbNH, "'E" and which cor- 
relate with the NH(A:3S-) + Rh(^S) dissociation limit. 
At nonlinear geometries (with Cg symmetry), these two 
states become '^A" and '^A" . These states are princi- 
pally bound by dispersion forces. However, as mentioned 
above, there is also an ion-pair threshold Rb"'"(^S') + 
NH~(^H) less than 4 eV above the neutral threshold. 
The resulting ^H ion-pair state of RbNH is subject to 
the Renner- Teller effect and splits at nonlinear geome- 
tries into two states with the unpaired electron (hole) 
either in the triatomic plane (^A') or perpendicular to 
it i^A"). There are thus two "^A" states, one of ion-pair 
character and one dispersion-bound, which form a con- 
ical intersection at linear geometries, where one has ^H 
and the other symmetry. [This is a conical intersec- 
tion on the two-dimensional surface obtained when the 
NH bond length tnh is fixed; in three dimensions, the 
two surfaces intersect on a seam at linear geometries.] 

The influence of ion-pair states on collisions of alkali 
metal atoms and similar species has been called the "har- 
poon" mechanism If the electron transfer takes place 
at long range, as is often the case, large collision cross 
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sections result. The range of the harpooning is deter- 
mined by the position i?x of the crossing point (conical 
intersection). This can be estimated from 



AEo 



(1) 



where AEq is the difference between the Rb ionization 
energy and the NH~ vertical ionization energy (NH elec- 
tron affinity). The experimental values for the Rb ion- 
ization potential and the NH vertical electron affinity are 
4.177 eV [l^ and 0.387 eV respectively; the latter is ob- 
tained from the adiabatic NH electron affinity (0.381 eV) 
[T?! and a calculated value of the difference between the 
zero-point energies of NH and NH~ (0.006 eV) This 
gives an estimate of i?x ~ 3.79 A. Because in NH~(^H) 
the negative charge is concentrated on the nitrogen atom, 
Rx represents the distance between N and Rb. 

The nature of the dispersion-bound and ion-pair states 
depends on the electron configuration of the molecular 
partner. For Rb-OH, for example, the dispersion-bound 
states are ^H and ^H at linear geometries, splitting into 
^A", ^A' and '^A" at non-linear geometries, while the 
ion-pair state is a closed-shell singlet (-'^E+ or ^A'). As for 
Rb-NH, therefore, there is no ion-pair state for the high- 
est allowed spin multiplicity. For Rb-CH, the dispersion- 
bound states are the same but there are many more ion- 
pair states arising from a tt^ configuration: {^A"), 
(}A' and ^A") and {^A'). The only dispersion- 
bound state that does not have an ion-pair counterpart 
in this case is the ^A' state. 

The high values of the dipole polarizabilities for alkali- 
metal atoms imply large dispersion forces. Because of 
this, the crossing point is in the classically allowed re- 
gion of the dispersion-bound state for many molecule - 
alkali atom systems. This is the case for RbNH. The 
"harpooning" is thus a barrierless process and no tun- 
nelling is required for the charge transfer. 

To obtain a more quantitative picture of the RbNH 
electronic states, we have carried out ab initio calcula- 
tions at the CASSCF(10,3)/MRCI level (complete active 
space self-consistent field / multireference configuration 
interaction). For N and H, we used the aug-cc-pVTZ ba- 
sis sets of Dunning [l9| in uncontracted form, and for Rb 
we used the quasirelativistic small-core ECPIOMWB ef- 
fective core potential of Leininger et al. ^20] , augmented 
by the valence basis set from Ref. IIJ. The (10,3) active 
space covers all the occupied orbitals on N and H and 
all the large-valence and 5p orbitals on Rb. Only the 
electrons in the Is orbital on N were excluded from the 
correlation calculations. All the calculations were per- 
formed using the MOLPRO package fS]. 

Figure n shows the resulting curves for Rb-N-H angles 
of 180° and 170° (for fixed NH bond length). The mini- 
mum on the quartet surface is at a depth of 0.078 eV at 
''RbN = 3.26 A, while that on the doublet ion-pair sur- 
face is at a depth of 1.372 eV at rRbN = 2.53 A. At linear 
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FIG. 1: CASSCF(10,3)/MRCI potential energy curves of 
RbNH witli fixed tnh = 1.0308 A and fixed Rb-N-H angle of 
180° (a) and 170° (b). (a) At linear geometries the crossing 
point is at 3.77 A. (b) At non-linear geometries the crossing 
is avoided. 



geometries the ^H and ■^E^ curves cross at i?x = 3.77 A. 
At non-linear geometries the crossing is replaced by an 
avoided crossing. Note that the minimum of the 
curve is at a shorter Rb-N distance (around 3.4 A) than 
the crossing point, so that at 170° the l^A" curve has 
just one minimum. 

The collisional processes in Rb-NH mixtures will de- 
pend greatly on the magnetic states that are present. If 
both NH(A^E~) and Rb(^S') are in spin-stretched states 
(|A/i?| = F = Fmax), coUisions will take place on the '^A" 
surface. If the NH(A'^S~) molecules are initially in their 
ground rovibrational state (z; = 0, A'^ = 0, J = 1, M = 
1) 14], only elastic and reorientation (spin-changing) 
atom-molecule collisions can occur. Under these circum- 
stances, trap loss will be caused only by spin-changing 
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collisions and by three-body recombination processes. 

If the species present are not in spin-stretched states, 
collisions will also involve the doublet surfaces. In this 
case, the harpoon mechanism may be important and may 
enhance the rates of reorientation collisions. In addition, 
the transient collision complexes formed during the col- 
lision will be highly polar, and may be expected to have 
enhanced interactions with other species. Three-body 
collision rates may thus be large. Collisions with third 
bodies might also result in vibrational quenching and for- 
mation of bound RbNII(^n) molecules with ionic char- 
acter, and consequently in trap loss. However, because 
of the short-range nature of this type of three-body re- 
combination (which can occur only when the distance 
between the N and Rb is smaller than Rx), it cannot be 
described in terms of the fragment-fragment scattering 
lengths and a three-body parameter. 

On the other hand, the harpoon mechanism followed 
by, say, photoassociation offers a possible means for pro- 
duction of extremely polar molecules, with electric dipole 
moments well over 5 D (the value of the electric dipole 
moment of RbNII(2n) varies from 8.7 D to 15.3 D be- 
tween the potential minimum and the crossing point). 
Such molecules might provide a valuable probe of low- 
temperature physics of dipolar gases, because phenom- 
ena related to dipole-dipole interactions, such as linking 
of polar molecules |2^, will be significantly enhanced. 

The possibility of producing ultracold fermionic 
molecules and achieving a transition to a superfluid phase 
is particularly interesting. The molecule - alkali atom 
complexes have particularly large dipole moments, so 
at first sight are very promising candidates, but it is 
worthwhile to seek a more quantitative characterization 
of which dilute dipolar Fermi gases will be most suitable 
for transitions to the superfluid phase. 

Following Baranov et al. 8] , we assume that the mean 
dipole-dipole energy per particle is much smaller than 
the Fermi energy, 



(2) 



where n is the number density of the gas, d is the molec- 
ular electric dipole moment and e^? is the Fermi energy. 
Because the relation "much smaller than" is rather vague 
and difficult to work with, we introduce an upper bound 
ctTjCh for the ratio of nd^ to ep 
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< QfuCL < 1- 
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Since there appear to be severe technical limitations to 
the cooling of atomic Fermi gases below 0.2e p/kB 
(where ks is Boltzmann's constant), we introduce a lower 
bound for the ratio of the critical temperature for the 
BCS transition, Tc, to the Fermi temperature sp/kB, 



kBTc 

£p 



> Pbcs- 



(4) 



Baranov et al. 
for this ratio, 



have derived an analytical expression 
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where = ?i(67r^n)^/'^ is the Fermi momentum and aa = 
—2md^/{n^h'^) is the effective scattering length. 

We can combine the above expressions and derive lower 
and upper bounds for the product md'^n^^^, 



(6) 

For a given species (i.e, given m and d), Eq. © gives the 
condition for the number density at which the dipolar 
Fermi gas reaches the so-called ultracold limit, Eq. (jSl, 
and the transition temperature satisfies Eq. Obvi- 
ously densities closer to the upper bound would give a 
higher Fermi temperature (and consequently higher crit- 
ical temperature) and might be preferable from an ex- 
perimental point of view. If for given aucL we introduce 
maximal density rimax by the equation 



aucL, 



then the ratio from Eq. Q is reduced to 

= 1.44 exp < >. 



ep 



12aucL J 



(7) 



(8) 



Values of aycL between 0.1 and 0.2 give kBTc/e-p ratios 
between 0.10 and 0.39. 

Fermionic ^*NH has an electric dipole moment of 
d — 1.389 D [2^, which gives a rather large absolute 
value of the effective scattering length, = —876 A. 
Setting aucL = 0.1 leads to a "reasonable" density, 
nma,x = 5.4 X 10^^ cm~'^. Consequently, we obtain a 
very favorable Fermi temperature ep/ks = 757 nK and 
critical temperature Tc — 80 nK. 

In the case of **^Rbi''NH(2n), the large dipole mo- 
ment and large mass give an enormous absolute value of 
the effective scattering length, = —234487 A. Setting 
Q^UCL = 0.2 gives a very low density, n„iax = 2.3 x 10®, 
which combined with the large mass gives tiny values, 
6.2 pK and 2.4 pK, for the Fermi temperature and crit- 
ical temperature respectively. That means that, despite 
the very large dipole moment, it would be experimentally 
difficult to achieve the transition to a superfluid phase in 
a dilute gas of ^'^Rb^'^NH. 

The above discussion demonstrates that the combina- 
tion of large mass and large electric dipole moment is 
undesirable for production of superfluid dipolar Fermi 
gases. Large mass implies low Fermi temperature, which 
consequently makes the superfluid transition tempera- 
ture very low. Table ^ summarizes the situation for 
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TABLE I: Dipole moment d (D), effective scattering length (A), maximal density Wmax (cm Fermi temperature £p/kB 
(nK) and critical temperature Tc (nK) for selected fermionic molecules in ground electronic states. 





d 






"UCL = 0.1 






QuCL = 0.2 




"max 




Tc 


n 


£p 1 kg 


Tc 


^Li^^Na 


0.46" 


-186 


5.7x10" 


8.7x10^ 


917 


4.6x101^^ 


3.5x10* 


1.4x10^* 


"NH 


Lssg' 


-876 


5.4x10^^ 


757 


80 


4.3x10" 


3030 


1180 




l.S'^ 


-1432 


1.2x10^2 


202 


21 


9.9x10^^ 


809 


315 


Hi2c"N 


2.98'' 


-7258 


9.5x10^ 


6.1 


0.65 


7.6x10" 


25 


9.6 


^Li«K 


3.41" 


-16919 


7.5x10** 


0.65 


0.068 


6.0x10^* 


2.6 


1.0 


^Li^^NH 


5.r 


-17335 


7.0x10** 


1.32 


0.139 


5.6x10^* 


5.3 


2.1 




8.72^' 


-234487 


2.8x10^ 


0.0016 


0.0002 


2.3x10^ 


0.0062 


0.0024 



"Reference 2^ 
'Reference 2^ 
'^Reference 27 
"^Reference 26 

"^Present work; RHF value at RCCSD(T)/cc-pV5Z optimized lin- 
ear geometry r^H = 1.014 A and tljim = 1.695 A. 
•''Present work; CASSCF(10,3)/MRCI value, see text for details. 



several fermionic molecules that are of interest to the 
ultracold community. The maximum number densities, 
Fermi temperatures and critical temperatures are listed 
for two values of aucL, 0.1 and 0.2. Table illustrates 
that the dipole moment and molecular mass are both rel- 
evant characteristics when considering transitions to su- 
perfluid phases in dipolar Fermi gases. Lighter fermionic 
molecules with a reasonably large dipole moment, such 
as ^^NH or ^Li^^Na, are more suitable for transition-to- 
superfLuid experiments than heavy molecules with very 
large dipole moments. 
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